Preprint typeset using I^Tf7)X style emulateapj v. 8/13/10 



LAGRANGIAN COHERENT STRUCTURES IN NONLINEAR DYNAMOS 

E. L. Rempel 1,2 , A. C.-L. Chian 1 ' 3 ' 4 , and, and A. Brandenburg 5 ' 6 

ABSTRACT 

Turbulence and chaos play a fundamental role in stellar convective zones through the transport of 
particles, energy and momentum, and in fast dynamos, through the stretching, twisting and folding 
of magnetic flux tubes. A particularly revealing way to describe turbulent motions is through the 
analysis of Lagrangian coherent structures (LCS), which are material lines or surfaces that act as 
transport barriers in the fluid. We report the detection of Lagrangian coherent structures in helical 
MHD dynamo simulations with scale separation. In an ABC-flow, two dynamo regimes, a propagating 
coherent mean-field regime and an intermittent regime, are identified as the magnetic diffusivity is 
varied. The sharp contrast between the chaotic tangle of attracting and repelling LCS in both regimes 
permits a unique analysis of the impact of the magnetic field on the velocity field. In a second example, 
LCS reveal the link between the level of chaotic mixing of the velocity field and the saturation of a 
large-scale dynamo when the magnetic field exceeds the equipartition value. 

Subject headings: magnetohydrodynamics (MHD) — dynamo — chaos 



1. INTRODUCTION 

The equipartition-strength magnetic fields observed in 
planets and stars are the result of a dynamo process, 
whereby kinetic energy from the motion of a conducting 
fluid is converted in to magnetic energy (Brandenburg & 
Subramanian|2QQ5 ) . Initially, a weak seed magnetic field 
B undergoes a linear growth in the so-called kinematic 
dynamo phase until B is strong enough to impact the 
fluid velocity u. Eventually, the magnetic energy sat- 
urates due to nonlinear effects. In a fast dynamo, the 
growth rate is positive and non-vanishing even in the 
limit where the magnetic Reynolds number tends to in- 
finity. It is known that the growth of the magnetic energy 
in fast dynamos is related to the presence of Lagrangian 
chaos in the velocity field, i.e., scalar quantities passively 
advec ted by the flow (passive scalar s) exhibit chaotic mo - 



tions (Childress & Gilbert 1995 Balsara & Kim|2005) 



As B grows, it may suppress this chaos due to backre- 
action in the velocity field via the Lorentz force, le ading 
to the nonlinear saturation of the magn etic energy ( |Cat-| 
taneo et al. 1996 Zienicke et al. |1998 ). A comparison 
between the chaoticity of the velocity field during the 
growth an d saturation phases of the dynamo was per- 
formed by |Brandenburg et aT] ( |1995[ ) . In this letter we 
reveal how magnetic fields can affect the transport of pas- 
sive scalars through the formation of transport barriers 
in the velocity field. 
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When probing turbulent transport of passive scalars, 
either Eulerian or Lagrangian tools can be employed. In 
the Eulerian approach, for a given velocity field, one 
can solve an advection-diffusion equation for the pas- 
sive scalar concentration from which a turbulent diffu- 



sion coefficient can be computed (Vincent et al. 1996). 
Moreover, instantaneous snapshots ol tracer and velocity 
fields can be used to e xtract coherent structures su ch as 
eddies and filaments ( jlsern-Fontanet et al. 2004). Al- 
ternatively, in the Lagrangian approach the dynamics of 
fluids is studied by following the trajectories of a large 
number of fluid elements or tracer particles. The La- 
grangian description has been gaining increasing atten- 
tion in the past de cade, for example in the study of com - 
pressible plasmas (|Schamel||2004| [Padberg et al. 2007). 



It has been suggested that Lagrangian tools are more ap 
propriate to analyze tracer patterns than their Eulerian 
counterparts, since they do not rely solely on snapshots of 
the velocity field, bu t measure transport p roperties along 
particle trajectories ( d'Ovidio et al.|2009 ). We adopt the 
Lagrangian approach to distinguish the transport prop- 
erties of three-dimensional (3-D) numerical simulations 
of compressible magnetohydrodynamic (MHD) dynamos. 
We detect Lagrangian coherent structures (LCS), which 
are material lines and surfaces in the velocity field that 
act as barriers to particle transport and have been de- 
scribed as the Lag rangian building blocks of turbulence 
( |Mathur et al.|2007| ). There are two types of LCS formed 
by distinct groups of fluid particles, one of them attracts 
other fluid particles and the other one repels them. These 
barriers have been used to study turbulence and trans- 
port in fluids and p l asmas through nume rical simulations 



([Green et all2007l iPadberg et al. | [2007|), lab oratory ex- 
periments (Ivoth et al.||2002 Mathur et al.||2007|), and 



observat ional data in oceans and a wide range ol appli- 



Peacock & Dabiri 20111" 



cations ( Sandules cu et al. 2007 Olascoaga et al. 2008 



Two dynamo models with helical forcing and scale sep- 
aration are used, the ABC-flow ( jChildress fc Gilbert| 
1995) and an isotropic flow driven by a forc e correspond- 
ing t o plane waves with random phases (Brandenburg] 
2001). In the ABC-flow, two different dynamo regimes 
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are investigated, a regime characterized by a robust spa- 
tially coherent mean-field and a regime with intermit- 
tent switching between coherent and disordered mean- 
field states. Here, the LCS reveal that the topology of 
transport barriers in the velocity field suffers a dramatic 
change when the magnetic field undergoes the transition 
from coherent to intermittent dynamo. In the second 
dynamo model, the randomly forced flow, we focus on 
the problem of the nonlinear saturation of the magnetic 
energy. We note a large difference between the pat- 
terns of transport barriers in the kinematic and saturated 
regimes. 

An ultimate goal of this project is to interrelate re- 
sults from the Eulerian and Lagrangian approaches. For 
example, within the Eulerian approach it has been possi- 
ble to compute mean-field dynamo transport coefficients 
and their ma gnetic quenching behavio r using the test- 
field method ( jBrandenburg et al.||2QQ8| . The quenching 
is being interpreted in terms ol a competition between 
kinetic helicity that results in an a effect and current 
helicity that produces a magnetic a effect. In an inho- 
mogeneous system, the local current helicity distribution 
results from a balance accomplished by magnetic helic- 
ity fluxes (see |Brandenburg fc Subramanian| ( |2QQ5|) , for a 
review). If there is a direct connection between LCS and 
the suppression of turbulent transport and, in particu- 
lar, the a effect, one might expect a certain correlation 
in the spatial patterns of these quantities. As a prepara- 
tory first step, we establish here the basic technique in 
the case of a homogeneous turbulent dynamo. 

2. INTERMITTENT DYNAMO IN THE ABC-FLOW 

The dynamo model adopted consists of the com- 
pressible MHD equations for an isothermal g as, as de- 
scribed by |Rempel et al. (2009a). The equations are 
solved with the PENCIL CODE Q in a box with sides 
L = 2tt and periodic boundary conditions, so the 
smallest wavenumber is k\ — 1. The sound speed is 
c s = 1, so our time unit is (c s &i) -1 and the unit of 
viscosity v and magnetic diffusivity r] is c s /k\. We 
add to the momentum equation an external forcing 
given by the ABC (Arnold-Beltrami-Childress) func- 
tion, f(x) = Af /y/3[(sm kfZ + cos kfy)x, (sin kfX + 
cos kfz)y, (sinkfy + cos kfx)z], where Af is the ampli- 
tude and kf the wavenumber of the forcing function. 
We use kf = 5 to obtain a separation between the 
energy injection scale and the scale of the box, and 
Af = 0.1, which ensures a root mean square veloc- 
i ty lirm s = (u 2 ) 1//2 < 0.4. Following |Rempel et al.| 
(2009a), a numerical resolution of 64 3 mesh points is 



chosen. The kinetic (Re) and magnetic (Rm) Reynolds 
numbers are based on the forcing scale, Re = u rms /vkf 
and Rm = u rms /r]kf, where v is the average kinematic 
viscosity and 77 the constant magnetic diffusivity. 

We fix v = 0.005, which in the absence of mag- 
netic fields corresponds to a weakly turbulent flow with 
Re « 16. For large values of 77, the seed magnetic field 
decays rapidly and there is no dynamo. After the onset of 
dynamo action at 77 « 0.053 (Rm « 1.5), the magnetic 
energy starts to grow at the expense of kinetic energy, 
until it saturates. Examples of magnetic structures are 



depicted in Fig. [I] for two values of 77 and different times. 
For r? = 0.01 [Fig. [jja)] there is a coherent large-scale 
By component that propagates along the z direction. For 
77 = 0.05 [Fig. [jjb)], the magnetic field displays an inter- 
mittent switching between ordered (t = 6500) and disor- 
dered (t — 9000) large-scale structures. The scale-bars 
reveal that B y at 77 = 0.01 is one order of magnitude 
stronger than at 77 = 0.05. 

A better understanding of the spatiotemporal dynam- 
ics can be obtained by computing B y , the ^-averages 
of B. The upper panel in Figure [2] shows the space-time 
evolution of B y and the mid panel shows the time se- 
ries of By at the point z = 0. The lower panel shows 
the spectral entropy S m (t), which is a measure of spa- 
tial com plexity computed from the power spectra of B v 
following |Rempel et aT] ( |2009b| and |Chian et al.| fl2010D . 
The left column refers to 77 = 0.01 and the right column 
to 77 = 0.05. The left column shows that the mean-field 
for 77 = 0.01 propagates like a robust spatially coherent 
dynamo wave. The corresponding velocity field displays 
a mean flow with propagating oscillations. The direc- 
tion of propagation is arbitrary and depends on the ini- 
tial condition, which shows that there is multistability 
in the system. For 77 = 0.05 (right panel) the mean- 
field is more fragile and there is on-off intermittency, 
with phases of spatially disordered patterns interspersed 
with phases of spatially coherent structures. We call the 
regime at 77 = 0.01 "dynamo wave" and at 77 = 0.05 "in- 
termittent dynamo". The time-averaged values of the 
spectral entropy are (S m (t)) t ~ 0.045 for 77 = 0.01 and 
(S m (t)) t « 0.33 for 77 = 0.05. 

The effect of the magnetic field on the velocity field 
and its transport properties can be studied using the 
maximum finite-time Lyapunov exponent (FTLE). The 
maximum FTLE, <7^ 0+r (cc), gives the finite-time average 
of the maximum rate of divergence or stretching between 
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Fig. 1. — Intensity plots of B y . (a) The propagation of a large- 
scale coherent pattern along the z direction for 77 = 0.01. (b) 
Switching between ordered (t = 6500) and disordered (t = 9000) 
patterns for 77 = 0.05. 
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Fig. 2. — (a) Space-time evolution of B y (upper panel), the time series of B y at z = (mid panel) and the spectral entropy Sm(t) 
(lower panel) for 77 = 0.01. (b) Same as (a) but for 77 = 0.05, displaying an intermittent dynamo. 



the trajectories of a fiducial particle at x(t) a nd its neigh- 
boring p articles from time t = to to t = to+r ( |Shadden et| 
al. 2QQ5| ). A positive o\ is the signature of chaotic stream- 
lines in the velocity field. Finite-time Lyapunov expo- 
nents are able to detect Lagrangian coherent structures, 
which are the time-dependent analogous of stable and 
unstable manifolds of invariant sets in time-independent 
velocity fields. For a three-dimensional time-dependent 
velocity field, regions of maximum material stretching 
generate local maximizing curves (ridges) in the FTLE 
field. Thus, repelling LCS (finite-time stable mani- 
folds) produce ridges in the maximum FTLE field in 
the forward-time dynamics and attracting LCS (finite- 
time unstable manifolds) produce ridges in backward- 
time flHaller||200l| |Shadden etHl|20Q5| |Padberg etH] 
2007]). 

Since backward-time integration of dissipative s ystems 



is a major problem due to numerical instabilities (Celani 
et al. 120041), we have to resort to interpolation of recorded 
data. A run from to — r to to ± r is conducted and 
full three-dimensional snapshots of the velocity field are 
saved at each dt = 0.5 time interval. Linear interpolation 
in time and third-order Hermite interpolation in space 
are used to obtain the continuous vector fields necessary 
to obtain the particle trajectories. For backward-time, 
the interpolated snapshots from to to to — r are used 
and the particle trajectories are computed with a fourth- 
order Runge-Kutta method. The choice of the spatial in- 
terpolation scheme may affect the local dynamics of indi- 
vidual particles, which can result in minor changes in the 
delimitation of some of the material lines detected. Here, 
we adopted a third order interpolation, which is the stan- 



FTLE (see Haller & Yuan ( 


2000); Shadden et al. (12005b; 


Padberg et aT (|2007|); Mat! 


kur et al. (|2007p; Mendoza & 



Mancho ( 2010[ )). We compared the results obtained for 
the ABC flow with the known results from the literature 
and they show excellent agreement. 

Figure [3] illustrates the difference between the LCS of 
the wave and intermittent dynamos. Figure [3^ a) is a vi- 
sualization of the z) components of the velocity field 
for r] = 0.01 at x = using the technique of line integral 
convolution, which shows the integral curves of (u y ,u z ) 
in different tones of gray. This snapshot was computed 
for to = 2000, when the magnetic energy of the dynamo 
wave has already saturated (see Figure 2(a)). The cor- 
responding LCS are shown in Fig. [3^b), where green 
and red lines represent the repelling and attracting ma- 
terial lines, respectively. The LCS represent the cr^ 0+r 
field computed with r = ±10. Figure [3^b) is plotted as 
a three-vector RGB image using Octave's imshow rou- 
tine, where the forward-time g\ field is stored in the 
"green vector" and the backward-time <j\ field in the 
"red vector" . Notice that the intersections between high- 
intensity red and green lines may produce yellow points. 
We stress that the important feature of these plots is not 
the absolute value of ai, but the ridges in its field, so 
the colormaps are normalized by the largest value of g\. 
Figures |3jc)— (d) plot the velocity field and LCS, respec- 
tively, for the intermittent regime at 77 = 0.05 (to = 2000 
and r = ±10). The LCS distinguish the dynamo wave 
and intermittent dynamo regimes quite well. This be- 
comes clearer in Fig. [4j which depicts enlargements of the 
rectangular regions in Fig. |3j For the dynamo wave [Figs. 
Qa)-(b)l, a large eddy in Fig. |4ja) is seen in the LCS plot 
01 Fig. H[b) as a "smooth" region with low level of parti- 
cle dispersion bordered by attracting and repelling mate- 
rial lines. The entanglement of attracting and repelling 
LCS is responsible for the transport of particles between 
eddies (in two-dimensional flows this transport mech- 
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anism is called lobe dynamics (|Rom-kedar 



Wiggins 

1990| ). The X-point marked with an arrow in Fig. TSb) 
specifies the location of a hyperbolic trajectory nearby a 
point where the velocity field is instantaneously zero in 
the (y, zVprojection [see the velocity field near the arrow 
in Fig. Ha)]. The material lines are cross sections of 
material surfaces, and trajectories approach the X-point 
along the green line and are repelled from it along the red 
line. For the intermittent dynamo [Figs. EJc)-(d)], the 
entanglement is much more complex, even tnough the ed- 
dies in Figs. [4ja) and[4^c) look similar. The Lagrangian 
plot unveils an intricacy of local structures which is not 
seen in the Eulerian frame. The arrows in Figs. [4^c)-(d) 
mark the location of an X-point whose time-dependent 
manifolds, or LCS, fill the phase space in such a way 
that a border of the eddy cannot be identified. In the 
intermittent regime, the numerous crossings between at- 
tracting and repelling material lines enhance transport 
between regions. This transport can be quantified by 
the maximum finite-time Lyapunov exponent <j^ 0+t of 
the particle trajectories. The mean value for the dy- 
namo wave at r] = 0.01, obtained from a PDF computed 
with the trajectories of 64 3 particles evenly distributed 
in the box, is <t^ 0+t ~ 0.25. For the intermittent dynamo 
at r] = 0.05, cr^ 0+r « 0.32 and, therefore, the resulting 
chaotic mixing is more efficient. Although the LCSs vary 
according to to, in general the LCS fields computed at 
r] = 0.05 display higher degree of complexity than at 
77 = 0.01. 

The enhancement in the flow's chaoticity when the 
magnetic diffusivity is increased from 77 = 0.01 to 77 = 
0.05 is the result of a reduction of the effect of the Lorentz 
force upon the velocity field. For 77 = 0.01, B rms « 0.27 
and w rms w 0.29, so the magnetic field can become strong 
enough to suppress Lagrangian chaos in the velocity 
field, inhibiting particle transport. As a result of lower 
chaoticity, as well as the backward transfer of magnetic 
energy from small to large scales due to kinetic helic- 
ity (a-effect) present in this system, the magnetic field 
B saturates in an ordered state with the scale of the 
box and the mean-field dynamics resembles a spatially- 
coherent propagating wave [Fig. [2|a)]. For 77 = 0.05, 
stronger magnetic diffusivity causes the magnetic field 
to be damped, with B rms « 0.076 and u rms « 0.38. A 
weaker magnetic field has small impact on the veloc- 
ity field and Lagrangian chaos becomes stronger, with 
enhanced particle transport and chaotic mixing. The 
chaotic motions of the flow carry the magnetic field lines 
and generate the disordered B field states. On the other 
hand, stronger chaos in the velocity field leads to en- 
hanced stretching, twisting and folding of magnetic field 
lines, which tends to cause the magnetic energy to grow 
QChildress fc Gilbert 1 1995| . The growth of B backreacts 
on the velocity field, suppressing chaoticity again and 
leading to the intermittent occurrence of ordered B field 
patterns observed in the intermittent dynamo of Figs, 
gb) andgb). 

3. NONLINEAR DYNAMO SATURATION IN THE B2001 
FLOW 

For thi s section we have p erformed computations de- 



scribed in |Brandenburg| ( |2001[ ) (hereafter, B2001), where 
the MHD equations are solved with a helical forcing with 
a time-dependent wavevector. At each time step, there 



is a random choice of wavevector with wavenumber kf 
around 5. The resulting flow is essentially the proto- 
type of the a 2 dynamo of mean-field dynamo theory. 
We have adopted run 3 of B2001, where v = 77 = 0.002 
and the Reynolds number based on the box size is about 
600 (Re = Rm = 18 as defined in Section 2 above) for 
a numerical resolution of 128 3 . Due to an inverse cas- 
cade of magnetic helicity discussed in B2001, the mag- 
netic field develops a robust spatially coherent mean-field 
pattern similar to the case 77 = 0.01 above [Fig. IlTa)]. 
Figure [5| a) shows the time-evolution of u rms and £? rms 
in log-linear scale, revealing the exponential growth of 
B rms in the kinematic phase, before saturation. Fig- 
ures [5| a), (b) show the LCS computed for r = ±10 at 
to = 100 (kinematic dynamo) and to = 1700 (satu- 
rated nonlinear dynamo), respectively. In the kinematic 
regime the patterns of material lines in the LCS plot are 
highly complex, and the chaotic tangle permeates the 
phase space, which favors the growth of magnetic energy. 
In the nonlinear regime, since B rms becomes consider- 
ably higher than u rms (superequipartition), the chaotic- 
ity of the velo city field is strongly decreased due to th e 
Lorentz force ( |Cattaneo et al.|1996||Zienicke et al.|1 998). 
The crossings between the main attracting and repelling 
lines are scarce, so there is comparatively little disper- 
sion of passive scalars and transport is inhibited. The 
level of chaotic mixing quantified by the average FTLE is 



for the saturated dynamo. 



0.34 for the kinematic dynamo and cr{ 



to+r 



0.18 



4. CONCLUSIONS 



We have shown that Lagrangian coherent structures 
(LCS) can be used for an in-depth exploration of par- 
ticle transport in 3-D MHD dynamo simulations. Our 



res ults ag ree with the previou s results by |Cattaneo et 
al. (1996) and Zienicke et al. (1998), who showed that 
the modification of the velocity field due to stronger B 
becomes clearer by examining the Lagrangian proper- 
ties of the flow as measured by the finite-time Lyapunov 
exponents (FTLE). Here, in addition to computing the 
forward-time FTLE field, the detection of attracting ma- 
terial lines as ridges in the backward-time FTLE field 
provides the pathways that are more likely to be followed 
by passive scalars in the fluid. Moreover, the superposed 
plots of both attracting and repelling LCS permit the 
identification of the principal mixing zones of the fluid. 

The two dynamo models adopted in this work exhibit 
weak turbulence, with reasonably low Reynolds numbers. 
Our goal was not to present state-of-the-art numerical 
simulations, but to introduce the LCS technique in the 
context of space/astrophysical plasmas using two impor- 
tant topics in the theory of nonlinear dynamos, namely, 
the onset of intermittency and the nonlinear saturation 
of the magnetic energy. Regarding the first topic, the 
connection between the Lorentz force and on -off dynamo 
intermi t tency in ABC flows was studied by Alexakis 



Pontis 



(120081); on-off intermittency has also been ob- 
served in laboratory experiments with a dynamo gen- 
erated by a flow of li quid sodium QRavelet et al.|[2008 1 



Monchaux et al.| [2009); besides, intermittent chaotic dy- 



namos have been suggested as the cause of the long pe- 
riods of low sola r activity in the solar cycle, known as 
grand minima (Spiegel 2009). Concerning the second 
topic, it is one of the fundamental questions in dynamo 
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Fig. 3. — (a) Line integral convolution plot of the (y, z) components of a snapshot of the velocity field at x = for the wave dynamo at 
r\ = 0.01. (b) The corresponding repelling (green) and attracting (red) LCS represented by material lines; (c) and (d): same as (a) and 
(b), but for the intermittent dynamo at 77 = 0.05. 

surface flows can also be inferred from helioseismic data 



theory and there is an extended list of papers that discuss 
the nonlinear saturation of B in dynamo simulatio ns in 
periodic boxes with moderate Reynolds numbers (Bran- 
denburg et al. 1995; Catt aneo et a L 1996; Zien icke et al. 
19981 |Brandenburg||2001[ | Bra ndenburg Subramanian 



2005 



2009 



Kapyla & Brandenburg 2009; Cattaneo & Tobias 



The LCS method can be readily employed in a number 
of problems related to the turbulent transport of pas- 
sive scalars, including observational data, provided an 
estimation of the velocity vector field is available. Such 
estimations can be obtained from digital images using 
tec hniques such as the opti c al flow algorithm, employed 
by |Colaninno fc Vourlidas ( 2006 ) to extract the veloc- 
ity neT^nxonTT^^ ejections obtained 
with the SOHO LASCO C2 coronagraph. Solar sub- 



(Woodard 2002), thus LCS can aid the tracing of particle 
transport by turbulence in stellar interiors. 

In conclusion, a proper understanding of Lagrangian 
chaotic mixing is crucial for understanding the dynamics 
of nonlinear dynamos as well as the elaboration of models 
of stellar interiors that can correctly account for the tur- 
bulent transport of particles, energy and momentum in 
convective zones, and we believe LCS are an innovative 
tool that should be further explored in Astrophysics. 
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The arrows point to t = 100 and t = 1700; (b) LCS for the kinematic dynamo at t = 100; (c) LCS for the saturated dynamo at t = 1700. 
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